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<n : 

^ ' A spin-isospin dependent Three-Dimensional formalism based on the momentum vec- 

tors for the four-nucleon bound state is presented. The four-nucleon Yakubovsky equations 
with two- and three-nucleon interactions are formulated as a function of the vector Jacobi 
momenta. Our formalism, according to the number of spin-isospin states that one takes 

( ^^ ) into account, leads to only a strictly finite number of the coupled three dimensional integral 

equations to be solved. The evaluation of the transition and permutation operators as well 
as the coordinate transformations due to considering the continuous angle variables instead 

' _h ' of the discrete angular momentum quantum numbers are less complicated in comparison 

with partial wave representation. With respect to partial wave the present formalism with 
the smaller number of equations leads to higher dimensionality of the integral equations. 

■ We have concluded that three dimensional formalism is less cumbersome for considering the 
' three-nucleon forces. 

lE' 

§1. Introduction 

(N ■ 
> ' 

The evaluation of the four-nucleon bound state properties because of the presence 

■ of the fourth nucleon is a challenging task both from the formalism aspect as well as 
the numerical one 1 )-. 19 ) The binding energy of the a-particle with both 2 + 2 and 
3 + 1 structures needs both realistic two- and three-nucleon forces in order to get 

. closer to the experimental number. 14 ) 

In order to introduce the formalism of four-nucleon bound state we choose the 
three dimensional (3D) approach. Let us seek the answer to this main question which 
^ . is our motivation to select this approach. Why do we use 3D instead of partial wave 

(PW) approach? Few-body calculations are traditionally carried out by solving the 
relevant equations in a PW basis. After truncation they lead to coupled equations on 
angular momentum quantum numbers. A few PWs often provide qualitative insight, 
but modern calculations need many different spin, isospin and angular momentum 
combinations. On the other hand the 3D approach replaces the discrete angular 
momentum quantum numbers with continuous angle variables and consequently it 
considers automatically all PWs. It may be useful to mention that in PW approach 
we should sum all PWs to infinite order, but in practice we truncate the sum to 
a finite angular momentum number which is dependent to the energy that we are 
working. It means that in higher energies we will need more PW components to 
obtain a convergence, whereas in 3D approach continues angle variables sum all PW 
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components to infinite order. So the number of equations in 3D representation is 
energy independent, whereas in PW representation it is energy dependent. It appears 
therefore natural to avoid a PW representation completely and work directly with 
vector variables. 

The motivation for developing this non PW approach is introducing a direct 
solution of the integral equations avoiding the very involved angular momentum 
algebra occurring for the permutations, transformations and especially for the three- 
body forces. 

So in contrast to the truncated PW approach, the number of equations in the non 
truncated 3D representation is energy independent. Therefore this non PW method 
is more efficient and applicable to the three- and four-nucleon scattering problems 
which consider higher energies than the corresponding bound state problems. To 
show this efficiency the leading order 20 ) and full formulation 21 ) of 3N scattering with 
the inclusion of 3NFs has been done. Certainly the 4N scattering formulation and 
calculations are major additional tasks to be carried out. 

We have recently applied the spin-isospin dependent 3D approach based on mo- 
mentum vectors to the 3N bound state. The Faddeev equations with NN interactions 
are formulated successfully as a function of the vector Jacobi momenta, specifically 
the magnitudes of the momenta and the angle between them, as well as the spin- 
isospin quantum numbers. 22 )' 23 ) This novel formalism, according to the number of 
spin-isospin states that one takes into account, leads to only strictly finite num- 
ber of coupled three dimensional integral equations to be solved. We have shown 
that this formalism for both 3 H and 3 He bound states yields the same number of 
coupled equations which for fully charge dependent case leads to only 24 coupled 
equations. As an application the spin-isospin dependent three dimensional Faddeev 
integral equations are solved with Bonn-B potential. 24 )' 25 ) Our result for 3 H binding 
energy with the value of —8.152 MeV is in good agreement with the achievements of 
the other partial wave based methods. The calculation of 3 H binding energy with 
most modern NN potentials, i.e. AV18 and chiral potentials, is currently underway. 
It should be mentioned that the calculation of 3 He binding energy, because of the 
presence of coulomb interaction, is more complicated in comparison to 3 H case and 
with the recent achievement considering coulomb interaction it can be included in 
the calculation more easily in future. 

We have developed in this article the 3D approach for the four-nucleon bound 
state studies with the realistic interactions. The present work lies a formal ground 
for numerical investigations which are planned. This will helps us to reach to the full 
solution of the four-nucleon bound state in a straightforward manner. We propose the 
spin and isospin degrees of freedom to be implemented into the three dimensional 
four-nucleon Yakubovsky equations. We work directly with the vector variables 
in the Yakubovsky scheme in the momentum space, which leads to two coupled 
sets of a very limited number of equations in the three vector variables for the 
amplitudes which greatly simplifies the calculations without using a partial wave 
(PW) decomposition. 

In the present article we also show that how three-nucleon forces (3NFs), which 
we have chosen for example the Tucson-Melbourne two-pion exchange 3NF for con- 
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sidering the full solution of the four-nucleon bound system, can be presented in very 
simpler formalism in comparison with the PW. 12 ** With regard to this simplifications 
the study of Tucson-Melburn 3NF effects in the 3 H binding energy is also achievable 
and currently is underway. 

This paper is organized as follows. In section [2] we briefly represent the cou- 
pled Yakubovsky equations for the four-nucleon bound state with two- and three- 
nucleon interactions. We also introduce the four-nucleon basis states in a realistic 
3D approach. In sections [3] and U] we derive the coupled Yakubovsky equations in 
the realistic 3D approach as a function of the momentum vectors as well as the 
spin-isospin quantum numbers with and without the 3NFs, respectively. We add ap- 
pendix|A]to show how would be less cumbersome in compare with the PW approach, 
if we evaluate the 3NF matrix elements for example for the Tucson-Melbourne 2n- 
exchange 3NF in the 3D approach. In section [5] we discuss about the number of 
coupled Yakubovsky equations in both the 3D and the PW approaches. Finally we 
summarize in section 

§2. The 4N Basis States in 3D Representation 

The bound state of the four nucleon (4N) system, in the presence of the three- 
nucleon forces, is described by the two coupled Yakubovsky equations: 

= G tP[(l - P 34 Mi) + |^ 2 >] + (1 + G t)G Q v£l\V), 
\fo)=G tP[(l-Pb)\ik) + \ih)], (2-1) 

where the Yakubovsky components |?/>i) and \ip2) belong to "3 + 1"(123, 4) and 

"2 + 2" (12, 34) partitions of the four particles respectively. 14 ) Here the free 4N 

propagator is given by Gq = (E — Ho) -1 , and Hq stands for the free Hamiltonian. 

The operator t is the NN transition matrix determined by a two-body Lippmann- 

Schwinger equation. P, P and -P34 are permutation operators. P = P12P23 + 

-P13-P23 permutes the particles in the three-body subsystem (123) and P = -P13P24 

(3) 

interchanges the two two-body subclusters (12) and (34). The quantity V^i defines 
a part of the 3NF in the cluster (123), which is symmetric under the exchange of 
the particles 1 and 2. This can be related by an interchange of three particles to the 

two other parts an d V^g that sum up to the total 3NF of particles 1,2 and 3: 

(1) (2) (3) (3) 
V 123 = V^i + V±;i + V" X 23 ■ For the 3NFs based on a meson-exchange picture, Vy 2 3 

describes the interaction induced by a pion interchanged between the particles 1, 2 

and on its path rescattered by the third particle, see Fig. [TJ Applying a combination 

of the transpositions to the set of the Yakubovsky components, one obtains the total 

wave function \Hf) as 

\W) = [1 - (1 + P)P 3 4](1 + PMi) + (1 + P)(l + P)\>h). (2-2) 

The antisymmetry property of the under the exchange of the particles 1,2 
and the 1^2} under separate exchanges of particles 1,2 and 3,4 guarantee that the 
is totally antisymmetric. According to the two types of partitions (123, 4) and 
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Fig. 1. 



Diagrammatic representation of the part of a 2-k exchange 3NF, like the TM force, and 



the definition of the momentum transfers Q and Q' within the two subsystems. 



(12,34), as shown in Fig. [21 we introduce two different sets of the 4N basis states 
in a 3D representation which are suitable to represent the Yakubovsky components 
I'i/'i) and 1^2) in the coupled equations f)2 • 1 1) : 

|u a) = |uiu 2 u 3 0:1234), 

I V (3) = I Vi V 2 V 3 /?1234 )- (2-3) 

The spin-isospin part of the basis states are: 

I , 1 S T \ 

I "1234 / = l«1234 a 'l234 /> 

= 1/^234 01234 >, (2-4) 

where 

I "1234 ) = I ((( s l S 2 )si2 S3>123 S A )S M S ) = I ((S12 ^123 ~)S M S ), 

\Pi23i) = I ((si s 2 )sn (S3 S4)S U )SM S ) = \ (s 12 s 34 )SM s ), (2-5) 

and the isospin parts of the basis states | 0^1234 ) and | /?i234 ) are similar to the spin 
parts. 

Each one of the basis states involves three standard Jacobi momenta Ui, 112 and 
U3 or vi, V2 and V3. 14 ) As indicated in the Fig. the angular dependence explicitly 
appears in the Jacobi vector variables, whereas in a standard PW representation the 
angular dependence leads to three orbital angular momentum quantum numbers for 
each kind of the basis states, i.e. £12,^3 and I4 or £12,^2 and Z34. 12 ) It indicates 
that in our 3D formalism there is not any coupling between the orbital angular 
momenta and corresponding spin quantum numbers. Therefore we couple the spin 
quantum numbers S12, S3 and S4 or S12 and S34 to the total spin 5 and its third 
component Ms, by using only one intermediate quantum number S123 for the first 
basis states: | ((S12 S3)si23 s^)SMs) or | (S12 334)5 Ms). For the isospin quantum 
numbers similar coupling schemes to the total isospin T Mt involve one intermediate 
quantum number £123: | ((£12 ^3)^123 t±)T Mt) or | (t\2 t^T Mt) ■ In order to be 
able to evaluate the transition and permutation operators we need the free 4N basis 
states I u 7 ) and | v 7 ) , where 

I 7 ) = I 71234 ) = I 7f234 7m4 ) , I 7f234 > = | m SJL m S2 m S3 m S4 ) , (2-6) 
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3D 




2+2: (12,34; 12) 




lui U 2 U 3 ((S12 S3)S 12 3S 4 )S M s ; ((t 1 2 t 3 )t r 23 t 4 )T M T > 

Ivi v 2 v 3 (si 2 s 34 )S M s ; (ti2 t 34 )T M T > 

PW 





lu i U 2 U 3 (((1 12 S 12 )j 12 (13 S 3 )j 3 )j 123 (U S 4 )j 4>J M j J ((t 12 t 3 )t j 23 t 4 )T M T > 

Ivi v 2 v 3 (((li 2 s 12 )ji 2 l 2 )j* 12 (l 34 s 34 )j 34 )JMj; (ti 2 t 34 )TM T > 

Fig. 2. Definition of the 4N basis states in both the 3D and the PW schemes which are constructed 
by 3 + 1 and 2 + 2 type of the Jacobi coordinates and the corresponding spin-isospin quantum 
numbers. 



where m Si ,i = 1,2,3,4 indicates the projection of the spin of each nucleon. The 
isospin part of the basis states | 7i23 4 ) is similar to the spin part. In this respect 
when we are going from one of the 4N basis states, | a ) and | ft }, to the free 4N basis 
states, I 7), we should calculate the following Clebsch-Gordan coefficients; 

(l\ a ) — 9*1 a = 571234 «1234 = 57x234 01234 ^71234 «1234 

= (m Sl m S2 m S3 m S4 \ ({s 12 ^)si23 ]r)SMs) 



1. 



1 



x ( m tl m t2 m t3 m u \ ((t 12 -)t 123 -) T M T ), 



(710: 



_ S T 

9*lP = 571234 /3l234 = #71234 /3l234 #71234 /3l234 

( m si m S2 m S3 m S4 \ (s 12 s 34 )S' M s ) 
x ( mti m t 2 m U m 4 1 (*12 hi)TM T ). 



(2-7) 



6 



S. Bayegan, M. R. Hadizadeh, and W. Glockle 



The introduced basis states are complete in the four-nucleon Hilbert space as: 
^ A |Ae>(A^| = l, ^ =Y,J D Z A = Y,J d 3 A, J d 3 A 2 J d 3 A 3 ,(2-8) 

where A indicates each one of the u and v vector sets and £ indicates a, (3 and 7 
quantum number sets. They are also normalized according to: 

(A£| A'O = <5 3 (A- A')5^>. (2-9) 

Clearly the basis states | u a } are adequate to expand the first Yakubovsky 
component and correspondingly the basis states |v (3) are adequate for the 
second one \ip2)- 

§3. Yakubovsky Equations in 3D Representation Without 3NFs 



Let us now represent the coupled equations (|2-ip . without the 3NFs, with respect 
to the basis states which have been introduced in Eq. (|2-3() : 



(ua|^i) =¥f {ua\G tP(l-P u )\u"a")(u"a"\ij 1 ) 

(v/3|V 2 )=yf (v/3|GotP(l-P 3 4)|u / a , )(u , a , |^ 1 } 
7a' 

+ T (v/3|G tP|v , /?'}(v'/?'|V2). (3-1) 



It is convenient to insert the free 4N completeness relations between the permu- 
tation operators, it results: 

(ua\^) = jF jT (ua|G tP|u'7 / )(u'7'|(l- J P34)|u ,/ a ,/ )(u"a"|V'i} 

V7' ifa" 

+ ]£" if (ua|G tP|u'7 , }(u'7 / |v'/3')(v / /3 / |V 2 ), 

< V /3|V> 2 )= ^ (v/3|Goi J P|v'7 , }(v'7'|(l-P 3 4)|u / a / )(u / a'jVi} 

ITy' Jo! 

+ JT {v(3\G tP\V f3')(V p'\ip 2 ). (3-2) 



For evaluating the above coupled equations, Eq. (|3-2p . we need to evaluate the 
following matrix elements: 

(ua\G tP\u' (3-3) 
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(v^lGotiW), (vP\G Q tP\v f 0'), (3-4) 

(u / 7'|l-P 3 4|u"a"), (3-5) 

(v , 7'|l-P 3 4|u'a'), (3-6) 

(u'YIv'/?'). (3-7) 

For evaluating the first term, Eq. (|3-3|) . we should insert again a free 4N com- 
pleteness relation between the NN i-matrix operator and the permutation operator 
P as: 

ru " 



u a |G tP|uV> = - g g M E^„ 

-° m 4m 3m 7"' ' 

x (a|7 //, )(u7 W |t|u' / 7 // )(u" 7 /, |i :, |u / 7'), (3-8) 



where the matrix elements of the NN t-matrix and the permutation operator P are 
evaluated separately as: 

( U 7 '" |i| u" 7" ) = <5 3 (u 2 - U 2 ) 5 3 (u 3 - U3) <5 m » 3TO « ^m' a ' 4 m« S m'^m'^ S m^m^ 



. . / W /// /// III \\ II II II II II \ 

4m 3m 



3v% 2ul , s 

E — — ^ — — ^, (3-9) 



I II II \r>\ J I \ xS/ II I \ x x 

(u 7 |P| u 7 ) = <5 (u 3 - u 3 ) d m ^ m ' s4 5 m » m j 4 



x I <5 3 « + ^ + jju 2 ) 5 3 (u' 2 ' - ui + 

x "m'^m' S2 "m'J 2 rn' S3 °m'J 3 m' 31 "m^m^ ""^Ki 

+5 3 K + iu' 1 -^u 2 )5 3 (u 2 ' + u' 1 + Iu 2 ) 



X ^ m " 1 m s 3 " l >2 m ! 1 ^ m s 3 m 's 2 rn t 1 m 't 3 ^ m " 2 m i 1 ^ m " 3 m t 2 



(3-10) 

For evaluation of the matrix elements of the permutation operator P we have 
used the relation between the Jacobi momenta in different three-body subsystems 
(312,4), (231,4) and (123,4). Inserting Eqs. $3^) and (|5TPjl into Eq. leads to: 



/ \n + v>\ 1 i\ <j3 ( U3 ~ u 3) x 
(ua\G tP\ U 7 ) = — — — ^ ^ W 6 m'» 4 rn' S4 K> t '> m > 4 

m 4m 3m 7 
X | 5 3 (u 2 - u'j + -U 2 ) ^m™m' Sl bm'llm,'^ 



x (ui m^m's" m£m£'|t(e)|— u 2 - u 2 m'^m^ m^m^) 
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+<5 3 (u 2 + u[ + ~u 2 ) S m m m ' s2 5 m '^ m ' t2 

w / HI III III lit I . t M i ' ' / / / \ I 

x (ui m Si m S2 m ti m t2 |t(e)|-u 2 + u 2 m S3 m sl m t3 m tl ) \. 

(3-11) 

Representation of the second terms, Eq. (|3-4p . follows the similar steps: 

(v p \G tP\v 0' ) = J2 < v igo^i v' y ) ( ) = E w < v \ G ° tp \ v ' v > , 

(3-12) 



where 



<v/?iGotivy) = — v2 1 v2 v2 y: i: 



E - 1L _ J3. _ ill — -ys>, 

m 2m m 7 

x ( 0|Y" ) ( v 7"' |t| v" 7" ) ( v" 7 " \P\ v' V >• (3-13) 

The matrix elements of the NN i-matrix and the permutation operator P are 
evaluated as: 

( v 7 W |t| v" 7" ) = <5 3 (v 2 - v 2 ) (5 3 (v 3 - Vg) ^ m » 3 m« ^m^m™ ^m^rn^ 

. . / /// III III III ' u / * \ | ff // // ff ff \ 

x (vim Sl m S2 m ti m t2 |t(e )| Vl m Si m S2 m 4l m 42 ), 
e* = P-^-^, (3-14) 

{ v" 7" |P| v' 7' ) = 5 3 (v'/ - v' 3 ) 5 3 (v 2 ' + v' 2 ) <5 3 (v 3 ' - v'O 

x ^m" m' "mi' m' ^m" mi ^m" m' ■ (3'15) 
Inserting Eqs. (EKLl) and (13^131) into Eq. <^W) leads to: 

5 3 (v 2 + v 2 )5 3 (v 3 -v' 1 ) 



(v^iGotPivy 



E _ vf _ j± _ 4 

m 2m m 



i Vl m" s [m"J 2 «|t(e*)|v , 3 m'^ ). (3-16) 



For evaluation of the third term, Eq. (|3-5p . we should use the relation between 
the Jacobi momenta in different chains (123,4) and (124,3), which results: 



( u' y |1 - P 34 | u" a" ) = J2 9Y>a« <5 3 (u'i ~ u?) 
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-PM - - j^) 5 3 (u' 3 - + iug) 5^ 243 



(3-17) 



And finally for the evaluation of the fourth and fifth terms, Eqs. (j3-6|) and ([3l 
we should use the relation between the Jacobi momenta in two naturally different 
chains (123,4) and (12,34). the result is: 

( u V I v' ) = gyp, 5 3 (v[ - ui) £ 3 (v' 2 + u' 2 + ^u 3 ) 5 3 (v 3 - l -W 2 + ^u' 3 )^348) 



( v' i |1 - P 34 | u' ol ) = 9T>a> ^ 3 (u'i " 



3^ 3^ A 2 T1234T1234 



/ 1 2,,/ ! 2^, ^ x 3/„/ , l, r / _ 

T1234T1243 f ' 



+<5 3 (u 2 + ^v 2 + ^v 3 ) 5 3 (u 3 + ^v 2 -v 3 )^ .» 



(3-19) 

By these considerations, in the following we evaluate the first and second Yakubovsky 
components separately. For the first component we can rewrite Eq. (|3-2p by consid- 
ering Eqs. (|3-17p and (|3- 18|) as: 

(ua|^i) = ^ U (ua|G tP|u' 7 ') ^ VW'^K-u'/) 

-* 3 (u 2 - iuS - jjug) * 3 (u 3 - u» + iu 3 ') ^ 2347 » 243 I 

x(u" a" 
fu ' /V 

+ X £ (ua|G tP|u'7') 

x 5 3 (v'i - ui) 5 3 (v 2 + u' 2 + ^u 3 ) ,5 3 (v 3 - ~u 2 + ^u 3 ) 
x(v'/3> 2 ), (3-20) 

by integrating over u" and v' vector sets in the first and the second terms respectively, 
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and considering Eq. (|3-lip . we obtain: 

T-f?' £ 3 (u 3 -u 3 ) ^ 



*4 *4 



m 4m 3m 7 



5 (u 2 - u[ + -U 2 ) 5 m w m ^ $rn'^m' ti 



. . / III III II II I . / \ I / / f / / \ 

x (ui m sl m S2 m tl m t2 |i(e)|— u 2 - u 2 m S2 m a3 m t2 m t3 ) 



+5 3 (u 2 + u[ + ~u 2 ) <5 m ^ m ' s2 ^ m j» m ; 2 



x (m m Si m S2 m tl m fe |t(e)|-u 2 + u 2 m S3 m Sl m H m tl ) 



9-y'a" (u'a" 

a" 

^ 18 1 

- <? 7 "a" ^ 2347 (' 243 « -u 2 + -u 3 u' 2 - -u' 3 a" IV'i) 



and by integrating over u' x and u 3 vectors 
(ua |V>i) = / d 3 u 2 ^ 9ay 



(3-21) 



8m'" m' 8 m '"m' 

s 4 s 4 £4 £4 



m 4m 3m 



$m'J' 3 m' n Sm'^m^ ("i m Sl m S2 m tl m t2 |t(e) | — u 2 

< \ Y 9 -y' a " ( U2 + ^ U2 u ' 2 Us a " 

I a" 

- ^ 07"a" 6 Y 12M Yi243 ^ U2 + ^ U 2 ^ U 2 + ^ U 3 U 2 - ^ U 3 «" 

a" ,7" 

1 2 12 1 

+ W ( u 2 + ~u' 2 - u 2 - -u 3 + -u' 2 - -u 3 \ip 2 ) J 



+ <W'm' <W'm' ( u l ™">" 2 KiX" W) b u 2 + u' 2 m'^m'^ m' ts m' tl 



I Y 9"f'a" ( -U2 - ^U 2 U 2 
I a" 



u 3 a" 



118 1 

Yj 9Y'a" fyvatTvt* ( ~ u 2 - o U 2 o U 2 + n U 3 u 2 ~ o U 3 «" l^l) 



2 z 3 
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i 2 12 

+ Yl ^'P' ( ~ U2 ~ 9 U 2 ~ U 2 - q U 3 + ~ u 2 ~ o U 3 $ \>h) 



3 2 



(3-22) 



In order to evaluate the Eq. (13-22H we consider the result of the permutation 
operators P\i and P34 action on the Yakubovsky components, the space and also the 
spin-isospin parts of the basis states as: 



Pl2^l) = 


-l^i), 




Pl2[V>2) = 


-1^2), 




^34^2} = 


-M, 




Pi 2 |u) = 


I -uiu 2 u 3 ), 




Pl 2 |v) = 


I -viv 2 v 3 ), 




P34IV) = 


|viv 2 -v 3 ), 




Pu\a) = 


f_\Sl+S2— 812 f Ul+*2-*12 1 a \ — I 


-) Sl2+il2 |a), 


Pn\P) = 


f\Sl+S2— S12 fXtl+t2—tM 1 0\ _ ^_ 


-) sl2+tl2 \P), 


Pu\P) = 


f\S3+S4— 534/ U3+*4-*34| 0\ — t 


-) S34+t34 |/3), 


P12 It) = 


m S2 m sl m S3 m S4 m t2 m tl m t3 m ti ) 


= 1 72134}, 


^34 1 7) = 


m sl m S2 m S4 m S3 m tl m t2 m ti m t3 ) 


= 71243)- 



(3-23) 

By taking the Eq. (|3 23|) into consideration the matrix elements of the NN t- 
matrix as well as the Yakubovsky components in the second term can be rewritten 
as: 

(ui m"[m™ m "[ m ZW)\\ n ^ + u 2 "4 "4 m 't 3 m 't 1 ) 

/ m m 111 _ /// 1 . / \ rt t-> 1 1 ' ' ' ' ' \ 

= (ui m Sl m S2 m tl m t2 \t{e)Pi 2 Pi2\-U2 + u 2 m S3 m Sl m t3 m tl ) 

I III III til til I . / \ O I / / / I I \ /n ni\ 

= {\iim Sl m S2 m ti m t2 \t(e)Pi2\ - -u 2 - u 2 m Sl m S3 m tl m t3 ), (3-24) 



X] 9-y'a" ( ~ u 2 ~ ^ u 2 u 2 u 3 a" 

a" 

= Y 9 7 'a" { - u 2 - ^u' 2 u' 2 u 3 a" \P12P12\1p1) 

a" 

= -(-) s " 2+ *i2 2 g ia ii { u 2 + ^u 2 u 2 u 3 a" |Vi), (3-25) 



118 1 

^ 5 7 "a» ^2347^243 ( ~ U2 - 9 U 2 ^ U 2 + g U 3 u 2 - ^ u 3 1^1/ 



9 7 "a" ^2347(' 2 43 ( ~ U 2 ~ o U 2 7T U 2 + | u 3 u' 2 - -U 3 a" |Pi 2 Pl 2 |^l 
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s 12+ r 12 



/I 9j"a" & 



T1234T1243 



118 1 
x(u 2 + -u' 2 -u 2 + -u 3 u 2 - -u 3 a" |Pi 2 Pi 2 |^i), 

9-y'P' { -u 2 - ~ u 2 - | u 3 + \u 2 - -U3 \ljj 2 ) 



(3-26) 



~ ( _U 2 - \ U> 2 ~ u 2 - | u 3 + ^U 2 - ^U 3 \P\2PvM2) 



1 ../ 2 



; / ^ 2 12 

)4 2 +*i2 ^ 5y/3 , ( U2 + _ u ' 2 _ u ' 2 _ _ U3 + _ u ' 2 _ _u 3 \P 12 P 12 \tfj 2 ) 



(3-27) 



After inserting Eqs. (ETM - dSTTl) in Eq. (g^g) we obtain: 



s 4 s 4 *4 M 



«f 2u§ 



m 4m 3m 

' 



(ui m"'?<" m"'r?"4"|i(e)|— u 2 - u' 2 m' S2 m' S3 m^rn^} <5 m /« mii ^"m^ 
: \ Y 9-r'a" ( "2 + ^u' 2 u' 2 u 3 a" 

I a" 

a", 7" 

+ X] W < U2 + ^ U 2 ~ U 2 - l U 3 + ^ u 2 - ? u 3 > 



+ (ui m sl m S2 m tl m t2 \t(e)P 12 \ - -u 2 - u 2 m' Sl m S3 m tl m t3 ) 6 n , n „,^ 



I a" 



, Ji \tii \ 1 
\s 12 +t 



^ 57" a" & 



7l2347l243 



12) (u 2 + -u 2 u' 2 u 3 a"|Vi> 



)^i'2+*'i'2 N ) ( U2 + \u' 2 -u' 2 + ^u 3 u' 2 - -u 3 a" 



2 z 3 



+ ^2 9i'P 
a 1 



>12+%2 I / u 2 + -u 2 - u 2 - -u 3 + -u 2 



1 , __, 2 



1 2 

2 -2 - 3 U 3 |V"2 



(3-28) 



The exchange of the labels m' Si ,m' ti to m' S2 ,m' t2 and reverse of it in the second 
term and consequently the following relations: 



g~f'a" — ► (— ) 



9-y'a" 
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1*10 



9-y'P* - (-) s ^+*i2 5y/3 , ( 3 . 29 ) 
lead to: 

f 3 f dm>J> 4 m' a4 $m' t "m' t 

in'" ^ m Am 3m 

X I X] g i a " { U2 + ^ U ' 2 U ' 2 U3 °" ^ 
I a" 



Y 9Y'a" ^ 2347 (' 243 ( "2 + ^U 2 -U 2 + -U3 U' 2 - ^U 3 «" | ^1 



2 2 3 2 9 

a", 7" 



(3-30) 



a.s: 



I 2 12 1 

+ Y 9j'p> { u 2 + ^ u 2 - u 2 - o u 3 + - o U 3 1^2) > , 
/3' ) 

For the second component we can rewrite Eq. (|3 2|) by considering Eq. ()3-19p 

(v/3i^) = ^; v (v/?iG tp|v' 7 ')^ vw^k-vo 

77 "TO 37 



5 3 (u' 2 + |v 2 - fv 3 ) 5 3 (u 3 + iv 2 + v 3 )5 



2 "2 I v 3;-7i2347l234 



obtain: 



+«5 3 (u 2 + 2 --v' 2 + ^v 3 ) 5 3 (u' 3 + iv 2 - v 3 ) ^ 2347 » 243 I X ( u' a' 

+ X (v/?|GotP|v , /3')(v / /3 , ^2). (3-31) 

In this stage by integrating over u' vector set and considering Eq. (|3-16j) we 
in: 

m 2m m 

X Y< 9 t 3 "f'" kjKj^mij ^m^ra;'_Jm; 2 
y 

. , / ill 111 111 111 1 . / * \ 1 / / / / / \ 
x {vxm sl m S2 m tl m t2 \t{e )|v 3 m S3 m S4 m t3 m u ) 

x j X] 9r r' a ' ( v 'i ~ ^ v 2 + ^ v 3 ~ ^ v 2 - v 3 «' IV>1> 

- Y 9Y'a' ^ 2347 » 43 ( vi - -v 2 - -v 3 - \w' 2 + v 3 a' 
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+J2 syW^hM \, (3-32) 

and in the next stage by integrating over and v 2 vectors we arrive to the following 
equation: 



hi mi 



3 / "'8 3 "*sx "'S4"'S2 "H 3 "Hi "H4, 7 'H 2 



7 >7 m 2m m 

x ( vi m"[m'" 2 m'l[w!" 2 \t(e*)\v' 3 m' S3 m' S4 m' t3 m' t4 ) 

{2 2 1 

Y 9~,>a> ( v 3 -v 2 + -v 3 -v 2 - v 3 a' 
a' 

~ 2 5 7; 243 «' < V 3 ^ V 2 - ^ v 3 + V 3 a' |V»1> 

a 1 

+ ^ 5 7'/3'( V 3 - V 2V 3 /3'|^2) >• (3-33) 

Therefore by considering the result of the permutation operators P i2 and P34 
action, Eq. (|3-23p . the matrix elements of the NN t-matrix as well as the Yakubovsky 
components, under the exchange of the labels m' S3 ,m' ta to m' S4 ,m' t4 and reverse of it 
with changing v 3 to — v 3 , can be obtained as: 

( Vl m™m% mZmZ\t(e*)\v' 3 m' S3 m' S4 m' t3 m[ 4 > 

— -> { vi m"[m" s ' 2 m£m£'|t(e*)| - v 3 m' S4 m' S3 m' u m' t3 ) 

1 111 111 in *\ d r> I ' ' ' ' ' \ 

= (vim Sl m S2 m tl m t2 |t(e )P i2 P i2 | - v 3 m S4 m S3 m u m t3 ) 

= (vim^m' s / >;>;; / |t(e*)P 12 |v 3 m; 3 m^m; 3 m; 4 ) (3-34) 

Y g ^' a ' ( V3 \ W2 + ^ V 3 ^ V2 - V 3 a ' l^l) 
a' 

2 2 1 

— » £ ^7 1243 «' ( v 3 3 V 2 - 3 V 3 2 V2 + v 3 a ' l^ 1 )' (3-35) 

a' 

Y H 2 43 a ' ( V3 ^ V2 ~ ^ V 3 ^ V2 + V 3 a ' l^l) 
2 2 1 

Y 9r r' a ' ^ V3 3 V2 + 3 V 3 2 V2 ~ v 3 a ' 1^' ( 3 ' 36 ) 

a' 

Y 9"fF ( V 3 ~ V2 V 3 Z 3 ' 1^2} » J] ^1243/3' ( V3 - V2 ~ V 3 Z 3 ' 1^2) 

= Y (-Y' 34+t ' M gyp* < v 3 - v 2 - v' 3 /?' 1^34^34^2) 
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= EH^s^ (-(-y'^A (v 3 -v 2 v 3 /3'k/> 2 ) 



= Sj'/3' (v 3 - v 2 v 3 /3' 1^2} • (3-37) 

With these considerations we can rewrite the Eq. ()3-33[) as: 



Vm^'m' Vm , ''m' Vm'J'mi ^m l J , m'+ 

w ^ — ^ — ^ 

7'i 7"' m 2m m 



. , / III III III III U I *\ T> I / ' ' ' ' \ 

x (vim Sl m S2 m tl m 42 |t(e )Pia|v 3 m S3 m S4 m t3 m t4 ) 



{ a' 



2 2 1 
v 3 g v 2 + -v 3 -v 2 - v 3 a' 



^234"' < V 3 I V 2 - I V 3 7J V 2 + V3 a' |-01> 



+ W (v 3 - v 2 v 3 /?' k/> 2 > i, (3-38) 



and finally considering the Eqs. (|3-33p and (|3-38p together leads to: 

1 f M , ^ S m%™^ <Wj< 2 $m!£m^ 



_./ _./// 



E _ ^ _ _ 4 

m 2m m 



x <vim^«|t(€*)(l -Pi 2 )|v 3 m^m' S4 m; 3 m; 4 > 

[ 2 2 1 

x < X] ^ V3 3 V2 + 3 V 3 2 V2 ~ Vs a ' 

I. «' 

2 2 1 

" f7( 2 34«' ( V 3 3 V 2 - 3 V 3 2 V2 + V 3 Q ' l^ 1 ) 



+ J2 9YP < v 3 " v 2 v 3 1^2} [ • (3-39) 

By introducing the physical representation of the NN i-matrix ( see appendix 
(B) of Ref. 23 )): 



a{pm Sl m S2 m tl m t2 \t{e)\ p' m' Sl m' S2 m' tl m' t2 ) a 

= (p m Sl m S2 m tl m t2 \t(e)(l - P u )\p' m' Si m' S2 m' tl m' t2 ) , (3-40) 



the final representation of the three dimensional Yakubovsky integral equations can 
be obtained by rewriting Eqs. (|3-30p and (|3 39|) : 
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X d 3 u' 2 &ry'" ^rn'l' A m' a Jm'i' 3 m' Sl ^ m 'l' A m[^m'^m' ti 

x a (ui m" s [m"J 2 m'l' 1 m'll\t{e)\—u 2 - u' 2 m' S2 m' s . } m' t2 m' t3 ) a 

X \ Y f7'a"( u 2 + ^ u 2 u 2 u 3 OL"\i>l) 
I a" 

" Y H***' ^ U2 + ^ U 2 ^ U 2 + ^ U 3 u 2 - ^ u 3a"|V>l> 
a" 

I 2 12 1 

+ Y ( U 2 + 2 U 2 ~ U 2 - o U 3 2 U 2 - 3 "3/^2) >, 

/3' J 



1 

yp\ih) = — ^-^-2 — ^ 

£J — _L _ _2_ _ 111 

m 2m m 



//,;'m' ^m'/'m' "m' t "m' t ^m'/'mL 



X / c*V 3 £ 

7' 7"' 

/ /// /// /// III I . / * M / / / / / \ 

x a {vim Sl m S2 m tl m t2 \t(e )|v 3 m S3 m S4 m t? m u ) a 
x I S 9 V<*' ( v 3 ^ v 2 + ^3 ^v 2 - Vga'lV'l) 

I. a' 

~ Y H243 a ' ^ V 3 ^ V2 ~ ^ V 3 ^ V 2 + v 3 «Vl 



Y 9 ~('P' ^ V3 ~ V2 ^Z 3 '!^) f- (3-41) 



It is worth to be mentioned that if we ignore the spin-isospin quantum num- 
bers we can easily reach the bosonic type of three dimensional Yakubovsky integral 
equations which are given in Ref. 18 ^ 

§4. Yakubovsky Equations in 3D Representation With 3NFs 

In this section we present the less cumbersome 3D representation of the Yakubovsky 
components with the 3NFs in compare with PW representation. The evaluation of 
the Yakubovsky components with the inclusion of 3NFs will be exactly the same as 
Eq. (|3-4ip except that an extra term with M^ 12 3 occurs i n the first component. This 
is 
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x ^(ua\vglm + (u a \tG v}£\V) j. (4-1) 

The matrix elements of the second term can be evaluated by inserting the suit- 
able completeness relations as: 

|-U "' 

(uq \tG vl~l\& x - ^ ^ 



EE? 

TT 37 ^a" 



^123 I 

x ( a | 7 ' ) ( u y |tGo I u w 7 " ) ( 7" I ) (u w a'" | vj% \ 9) 



EE^ 



07' a 57"a" 



# _ ill _ 3m 2" 2 _ 2u 'i' 2 



' ^ m 4m 3m 

,(3) 



x(u 7 'it|u w 7 ")(u'"a'"|^|^), (4-2) 



after evaluating the matrix elements of the NN i-matrix, Eq. ()3-9f) . and integrating 
over u'2 and u'^' vectors, we obtain: 

(»a|«w58i»>= E /"X" - *-»gr M 

7':7">a"' ^ m 4m 3m 

x ( ui m' sl m' S2 m' tl m' t2 \t(e)\u"' m'^m^m^m'Q 

x (u / 1 "u 2 u 3 a / "|F 1 ( ||<P'), (4-3) 

By using Eq. (|3-9p . the symmetry property of the 3NF and the anti-symmetry 
property of the total wave function under exchange of nucleons 1 and 2, one can 
write: 

(3)| l7 >\ _ [ rf.jrr Qi'oiWa." 



;ua\tG v£l\*)= [ d3 < 

_ J _ // . /// ** 



u'(' 2 3u| 2n| 
m 4m 3m 



x (uim Sl m S2 m tl m t2 |t(e)Pi2| - u x m S2 m Sl m t2 m tl ) 
x (_ui"u 2 U3« w (y/lllF), (4-4) 

under the exchange of the labels m" ,m" to m" 2 , m" 2 , reverse of it and changing 
, to-uf 



u'/' to -ui' we find: 



7' ,7", a'" " m 4m 3m 

I I I I I I if \T> \. Ill II II II II \ 

x (uim Si m S2 m ti m <2 | -t(e)Pi 2 |u 1 m Si m S2 m ti m t2 ) 

x «u 2 u 3 a'"|V^|!P). (4-5) 



3 
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Now we consider Eqs. (14-3P and (14-5P together to achieve: 

,3 /// fiy a 9j"a'" 



ua\tG v£ ) 3 \#) = ± Yl / dV i 



r ,r',a"' ' ^ m 4m ~~ 3m 

x (ui^ 1 m' S2 mJ 1 mJ 2 |t(e)(l-Pi 2 )|ui // m / s ' 1 m^m^m^) 

x « u 2 u 3 a"'\v}i\V), (4-6) 

by applying the introduction of the anti-symmetrized NN i-matrix, Eq. (!3-40p . we 
can rewrite the final representation of the three dimensional Yakubovsky integral 
equations, Eq. (|3 41|) . as: 

1 

(ua \rpi) 



u\ 3u\ 2u\ 
m 4m 3m 



J "is4 m s 4 ^rn'^m' n ^m^m^ ^m' t ^m' H 



/ d3u 2 Y 9aY" * 

x (ui m'^m^ m£'m£'|*(e)|-^u 2 - u' 2 m' S2 m' S3 m^m^ 

x | X] #7'o"( U 2 + ^ u 2 u 2 U 3 a"|V>l) 
I a" 

5, 7; 243 a" ( u 2 + 3 U 2 + g U 3 "2 - 3 U 3 a \W\ 

a" 

Y 2 12 

+ Yl W ^ U2 + 9 U 2 - U 2 - o U 3 ~ u 2 - o U 3 P'W, 



/3' 



+ (ua|t&V> 



1 „ f f 5 m' a3 m,'^m' aA m^J m [ 3 m'l^m[ i m'l A 



l',l",a.'" ^ m Am 3m 

x a (m m' Si m^ 2 m^m^ |t(e)|ui m'^m'^ m£m£) (ui u 2 u 3 a'" \V^iW) 



|T/(3)| 



(V/3|V2 



2m 



X / d 3 v' 3 Yj 5/37'" <5m™m' 3l <W^m' S2 ^m^m^ #m^? 
7', 7"' 

x o(vi m^m^ m1lm^\t(e*)\v' 3 m' S3 m' Si m' t3 m' u ) a 
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I a' 



2 2 1 
v 3 -v 2 + -v 3 -v 2 - a'lV'i) 



(4-7) 



To represent the generality of our 3D formalism we can simplify the Eq. (|4-7p 
to the bosonic case by switching off the spin-isospin quantum numbers, see Ref. 19 ) 
In order to show the efficiency of our formalism we have also chosen a realistic 2tt- 
exchange 3NF, i.e. TM, to evaluate the matrix elements of (ua \V\23 1^)- As a matter 
of reference we present this formalism in appendix (jA|) to indicate the simplicity of 
3D representation. 



§5. The Number of Coupled Yakubovsky Equations in both the PW 

and the 3D Formalisms 



In this section we discuss about the number of coupled Yakubovsky equations in 
both the 3D and the PW approaches. Already in the PW representation of the 3N 
bound state equations the number of channels must be high, e.g. of the order of 34, in 
order to achieve reasonably well enough converged energy eigenvalues. 24 ) In contrast 
to the 3N system the number of channels for the 4N bound state N = N a +Ng, where 
N a and Ng are the numbers of a and (3 quantum number combinations respectively, 
is in principle unlimited even if the 2N interaction is assumed to act only in a certain 
2N states. Nevertheless the approach to that unlimited number can be classified for 
instance in the following possibly useful manner. For fixed I4 and fixed total quantum 
numbers J and T the numbers of a-states is strictly finite once the 2N interaction is 
assumed to act only up to j]% ax ■ Correspondingly, the number of /3-states is strictly 
finite once I2 is fixed and again the 2N interaction is assumed to be zero in states for 
312 , J34 < jij ax ■ We display examples for those maximum N a and Np values in Table 
[5j Thus even if assuming only j'^ ax = 1, and restricting I2 and I4 to be at most 1, 
one leads to N = 60 channels, while the total isospin is restricted to be zero. 9 ) 

In Tables HT1 and Hill we present the number of the spin-isospin states for both kind 
of the Jacobi coordinates, a and (3, as well as the number of coupled Yakubovsky 
equations in our realistic 3D formalism. Since the angular momentum quantum 
numbers, i.e. l\2,h^4 an d ^i2)^2>^34> do not appear explicitly in our formalism, 
therefore the number of coupled equations which are fixed according to the spin- 
isospin states are strongly reduced. This is an indication that the present formalism 
automatically consider all partial waves without any truncation on the space part. 
Considering the spin-isospin degrees of freedom for 4 He one leads to 8, 38, 58, 18, 
87, 132, 22, 106 and 161 coupled equations for different combinations of total spin- 
isospin states (5 - T): (0 - 0), (1* - 0), (2* - 0), (0 - 1*), (1* - 1*), (2* - 1*), 
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Table I. The number of partial wave channels contributing to both kind of the 4N Jacobi coor- 
dinates for = + dependent on maximal values of intercluster orbital angular momenta U 
and total two-body angular momenta jij. For (a) and (b) parts the total isospin is restricted 
to T = and for (c) part T = 0,1,2. N is the total number of channels, where N a and Np 
are the number of channels corresponding to 3 + 1 and 2 + 2 partitions. The results of part (a) 
are given according to the notation of Ref., 9 ' where ji2, J34 < ji] ax and h,h < lT ax \ an d the 
results of parts (b) and (c) are given according to the notation of Ref., 14 ' where j'12, j'34 < jij ax , 
h,h,h < lT ax an d I12 + I3 + I4, I12 + '34 + h < It should be mentioned that for (c) part, 

N a = 4200 and Np = 2000. 



(a) Ref. 9) 












•max 
Jij 


N- 

J * (.1: 




N 


jmax rt 
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10 


10 


20 
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18 


18 


36 
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26 


26 


52 


jmax 












1 


34 


26 


60 




2 


66 


58 


124 




3 


98 


90 


188 


jmax 2 












1 


62 


34 


96 




2 


130 


98 


228 




3 


202 


170 


372 


jmax g 












1 


90 


34 


124 




2 


198 


122 


320 




3 


322 


242 


564 


jrnax ^ 












1 


118 


34 


152 




2 


266 


130 


396 




3 


446 


290 


736 


(b) Ref. 13 ' 










■max _ a 
Jij u 


jmax g 


imax -10 

L sum lz 


T = 


N = 1572 



(c) Ref. 14) 

j™ ax = 6 lT ax = 8 C = 14 T = 0, 1, 2 N = 6200 



(0 — 2*), (1* — 2*) and (2* — 2*) respectively. The star superscript indicates all spin 
or isospin states that we have taken into account up to a specific value. It is clear 
that in the 3D formalism for a fully charge dependent calculation there is only 161 
coupled equations, whereas in PW approach after truncation of the Hilbert space 
to T = there is 1572 coupled equation, see part (b) of Table and 6200 coupled 
equations for a fully charge dependent, see part (c) of Table [SJ So our 3D formalism 
leads to a very small number of coupled equations in comparison with the very large 
number of coupled equations in the truncated PW formalism. 

However, it should be mentioned that our formulation leads to coupled integral 
equations which depend on three vector variables for the amplitudes, whereas the 
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Table II. The number of the spin-isospin states for both kind of the 4N Jacobi coordinates in a 
realistic 3D formalism. N as (N aT ) and Np s (Np T ) are the number of the spin (isospin) states 
for 3 + 1 and 2 + 2 partitions correspondingly. 



(a) 3+1 partitions 



((512 |)S123 \)SM 3 
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5 = 0,1 


5 = 0,1,2 


((o pi s)o 
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((i P§ s)i 
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N as 2 11 16 
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((1 PI Pi 
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((1 PI Pi 
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N aT 
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(b) 2 + 2 partitions 








(sia s 34 )SM s 
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5 = 0,1 


5 = 0,1,2 


(0 0)0 
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1+0 


1+0+0 


(0 1)1 
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0+3+0 


(10)1 
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0+3+0 


(1 1)0 
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1+0 


1+0+0 


(1 1)2 





0+0 


0+0+5 




2 


8 


13 


(tl2 t34)T0 


T = 


T = 0,1 


T = 0,l,2 


(0 0)0 
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1+0 


1+0+0 


(0 1)1 





0+1 


0+1+0 


(10)1 





0+1 


0+1+0 


(1 1)0 
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1+0 


1+0+0 


(1 1)2 





0+0 
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PW formulation after truncation leads to a finite number of coupled equations in 
three variables for the amplitudes. 

The coupled Yakubovsky equations, Eq. (|4-7p . represent a set of three dimen- 
sional homogenous integral equations, which after discretization turns into a huge 
matrix eigenvalue equation. In order to solve these coupled integral equations di- 
rectly without employing the PW projections, one has to define a coordinate system. 
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Table III. The number of coupled Yakubovsky equations for a particle in a realistic 3D formalism 
according to the spin-isospin states (S — T) that we have taken into account. N — N a +N/3 is the 
total number of coupled Yakubovsky equations, where N a = N as x N aT and Np = Np s x N@ T 
are the number of 3 + 1 and 2 + 2 states correspondingly. The star superscript indicates all spin 
or isospin states that we have taken into account up to a specific value. 



(S-T) 


N a 


Np 


N 


(0-0) 


4 


4 


8 


(l*-0) 


22 


16 


38 


(2* - 0) 


32 


26 


58 


(o-r) 


10 


8 


18 


(r-r) 


55 


32 


87 


(2*-l*) 


80 


52 


132 


(0-2*) 


12 


10 


22 


(l*-2*) 


66 


40 


106 


{T _ T) 


96 


65 


161 



According to experience of 3 H binding energy calculations, Ref., 23 - 1 it is convenient 
to choose the spin polarization direction parallel to the z-axis and express the mo- 
mentum vectors in this coordinate system. As indicated in Ref. 19 ) generally one 
needs six independent variables, three magnitudes of the Jacobi momentum vectors 
and three angle variables, to uniquely specify the geometry of the three vectors. But 
in contrast to the bosonic case, Ref., 18 ) one does not have this freedom to choose 
the third vectors, i.e. 113 and V3, parallel to the z— axis and second vectors, i.e. U2 
and V2, in the x — z plane. So the angular dependencies of spin-isospin-dependent 
Yakubovsky components will be more complicated rather than bosonic ones. The 
dependence on the continuous momentum and angle variables should be replaced in 
the numerical treatment by a dependence on certain discrete values. Let the number 
of these discrete points be denoted by Nj aco u and N ang i e corresponding to momen- 
tum and angle variables, the dimension of the eigenvalue problem in both PW and 
3D approaches is: 

N PW = Nf acobl x {N™ + N™) 

N 3D = Nf acobl x Nl ngle x (Nl D + N$ D ) (5-1) 

It should be indicated that in the formulation of the 4N bound state we can 
consider the parity and the time reversal invariance to cut down the dimension of 
the Yakubovsky components, see appendix ([B]). So according to our experience in 
four-bosonic calculations and also by considering the parity and the time reversal 
invariance one needs in 4N realistic calculations thirty grid points for Jacobi momen- 
tum variables and ten grid points for angle variables. By these number of mesh grids 
the dimension of the eigenvalue problem for a fully charge dependent calculation, 
i.e. for total spin-isospin states (S* — T*) = (2* — 2*), will be: 

N PW = 30 3 x (4200 + 2000) ^ 1.7 x 10 s 

N 3D = 3Q 3 x 1Q 3 x (gg + g5 ) ^ 4.4 x 10 9 (5.2) 
Thus, the price for the smaller number of equations in 3D representation is the higher 
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dimensionality of the integral equations. In other words, algebraic simplification is 
achieved by a more involved numerical scheme. As indicated in ref 18 ' for bosonic case 
these dimensions reduce to the values N PW = 2.7 x 10 6 and iV 3D = 2.7 x 10 7 . As it is 
clear from the close similarity of dimension values and also based on the experience 
of simpler bosonic calculations, we expect that for more rigorous nucleonic case the 
numerical calculation could be achievable. 



§6. Summary and Outlook 



We propose a new representation of three dimensional Yakubovsky equations for 
the four-nucleon bound state including the spin and isospin degrees of freedom in the 
momentum space. This formalism stays closely to the bosonic structure where the 
spin and isospin degrees of freedom are ignored. This is an important step forward 
since the formulation based on partial wave decomposition, which includes the spin 
and isospin degrees of freedom, after truncation leads to two coupled sets of a finite 
number of coupled equations in three variables for the amplitudes. In contrast our 
3D formulation leads to two coupled sets of a strictly finite number of equations in 
three vector variables for the amplitudes. The comparison of 3D and PW formalisms 
shows that our 3D formalism avoids the very involved angular momentum algebra 
occurring for the permutations and transformations and it is more efficient and less 
cumbersome for considering the three-nucleon forces. This formalism enables us to 
handle realistic 2N and 3N potentials with all their complexity in four-nucleon bound 
state calculations. This work provides the necessary formalism for the calculation of 
four-nucleon binding energy which is under preparation. 
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Appendix A 

/(3)i 
123 I 



(3) 

The Evaluation of ( u a 

1^123 1 3 *) f° r the ™ 271 -exchange 3NF 



A.l. Preparation of the TM 3NF for the 3D Representation 

For the evaluation of the coupled Yakubovsky equations, Eq. (j4-7|) . matrix 
elements of the form (uo |^i23 1^) nee d to be calculated. In this section we evaluate 
these matrix elements in our 3D approach for the TM 27r-exchange 3NF. In the 
notation from Fig. HJ this force is given by 26 ) 

V 123 — V °2 • *4 01 * <4 



Q 2 + ml Q' 2 + ml 
x [n • r 2 [A + B Q • Q' + C (Q 2 + Q' 2 )) + D r 3 • n x r 2 a 3 ■ Q x Q'] , 

(A-l) 

with the momentum transfers Q and Q'. The cij's are Pauli spin matrices and the 
form factors are chosen to be F(Q 2 ) = ( ^ 2 ~gf ) 2 with the cut-off parameter A. One 
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distinguishes four terms in the TM force, the so called A-, B-, C- and L>-term. 
The constants A, B, C and D are determined by the low energy theorems. 27 ** The 
momentum transfers are given by 

Q = ki - 14 = { (+Ui - ~U 2 ) - (+u'i - \u' 2 ) 1 = {u 2 - U 2 }(231,4) 5 

I / / J (123,4) 

Q' = k' 2 - k 2 = { (-lli - 77 u 2) - (~ u l - 75 u 2) \ = {u' 2 - U 2 }(312,4) ,(A-2) 

I 1 £ J (123,4) 

as indicated in Fig. (TJ In order to be able to evaluate the TM 3NF matrix elements 
we should evaluate the scalar product of the spin and momentum transfer vectors, 
<j\ ■ Q and cr 2 • Q', the scalar product of both the momentum transfer vectors, Q • Q', 
the scalar triple product of the spin and momentum transfer vectors, 03 • Q x Q', 
as well as the isospin operators. Toward this aim we can rewrite the spin-space 
operators in the form that can be evaluated easily in our 3D representation as: 

o-i • Q = Q a\ ■ Q, 

a 2 -Q' = Q' a 2 • Q', 

Q Q' = QQ'i, 
03 • Q x Q' = |Q x q'| . qTq' 

= W a fe " 2^ (fT3 ' ^' )2 + Ya ** ■ ® UZ ■ ^ 

+ (as ■ ^ /)2 (fTs ' ^ )2 " ir a {as ' } ' (A ' 3) 



where 



7 = QQ' ,«= \/l-7 2 . (A-4) 
The scalar product of the spin-momentum vectors can be evaluated as 

<W>.Q|£<) =£^4> S (Q) 4: ms (Q) = 0i m ,, (A-5) 



where D^, mg is Wigner D-function which is defined generally as D s mt ms (q) = (zsm' s \qsm s 
The application of the TM 3NF to the total wave function \&) can be considered as 
sum of the the four independent terms: 



D D 



vSlm = £vf jw v®\9) = $>«>, (a-6) 



i=A 



where we have abbreviated the isospin operators by : 

Tx • r 2 i=A, B, C 



/« = T1 ' T2 (A-7) 

' T3 • T\ X T 2 1=D. V ' 
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A. 2. Explicit 3D Evaluation of the A-, B-, C- and D-Terms 

In the following we evaluate the matrix elements of ( u a terms. From Fig. Q] 

(3) 

we see that 1^23 

can be splitted into two parts and each part contains a meson 
exchange. The first meson is exchanged in the subsystem (31), where it is called 
for convenience subsystem 2, and the second is exchanged in the subsystem (23), is 
called 1. Since the structure of the 3BF is specified by two momentum transfers of 
consecutive meson exchanges, it is convenient to insert a complete set of states of the 

(3) 

type 2 between V^l and and another complete set of states of type 1 between 

(3) 

the two meson exchanges. Then the matrix element of V 12 3 is rewritten as 



Here the subscripts 1,2,3 of the bra and ket vectors stand for the different types 
of three-nucleon coordinate systems of (3 + l)-type fragmentation (ijk, 4). The 
coordinate transformation from the system of type 1 to one of type 3 can be evaluated 
explicitly as: 



The matrix elements of the isospin coordinate transformations can be rewritten 

as: 




(A-8) 



3 (ua | u'a')i = g, 



, 5 3 (u[ + i Ul + ^u 2 )5 3 (u 2 - Ul + iu 2 )5 3 K " u 3 ).(A-9) 




/// in \ i ii ii i in in \ i ii i t(€) i /// \ 

u a ) 2 = i(u a s \u a s ) 2 i( a T \B J \ a T ) 2 



(A-10) 



where 



!<u* ag | u'" a£ ) 2 = 9^111 5 3 « + + \^)5^ - u'/ + \^)5\^ - ug). 

(A41) 



The matrix elements of the coordinate transformations i ( a'^ | JW | a'^ ) 2 are given 
m Ref. 12 ) The isospin matrix elements have been derived in Ref. 29 ^ for the 3N system 
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and are given below for the 4N system for completeness: 



i(<4|/»K' 



St"t"'Sm^m^'Shu, C ( _ 6) (- 
1 



2 

t" 

t 23 




,J:i4/ *23*31 



i=A, B, C 



in 

23"31 



5 T » T '''5 M » M »'^» 3i ^» 2 24z (-) 



(A42) 



r 231 
1 
2 

,/// 

r 31 




1 I f // 

2 2 t 23 



i=D. 



In the following steps the matrix elements for the different V's are evaluated. 
In the first step we evaluate the matrix elements of V 3 \. Since both pion-exchange 
propagators in the 3NF term only depend on the momentum transfer in a two-body 
subsystem, as indicated in Eq. (|A-2p . and also because of the separation of the isospin 
parts of 3NF, the matrix elements of V 3 \ can be written as: 

1 (u'a' \V 3 f | u" a" h = 5 3 (u'i - u?) 5 3 (u' 3 - ug) S^i 1 ( u 2 a' s \V 3 f | ug a" s > lf 

(A43) 

and the spin-space parts can be more simplified for the A-,B-,C- and Z)-terms sep- 
arately as: 



/ i i I t AA.B^C) I // ii \ i i x x x 

1 \ u 2 a S I V 31 I u 2 a S ) 1 = 2^ Z> 9 a' s ~(' s 9 a'^'S ° m 's 2 < 2 d < 3 < 3 °™>'s 4 < 4 



x i(u 2 m si \V 31 |u 2 m si )i, 



i ( u 2 a' s | V 31 | u' 2 ' ag ) i = gl> s y s #i» 7 » <5< 2 m» 2 ^m' S4 mi' 4 

7^ 7^ 

/ / / / nr(D)\ ii ii ii \ 
x i(u 2 m Si m S3 \V 3 \ ; |u 2 m Sl m S3 )i, 

where the A-,B- and C-terms can be evaluated as: 



(A44) 



(A45) 



i(u 2 m' Sl |y 3 f |U2<)1= i 



' ' I 77V /02\ °"l " Q II II \ 

u 2 m si \F{Q ) Q2 + m2 Q\ u 2 m Si )j 



(u 2 



u 



/'I ^("2-0 
»>2 +m 2 I "2 a 2l W m , s m , L , 



Uo 



(A46) 



' ' \iAB)\ ii ii \ 
U 2«V IV31 '|Ua m 8l ) 



, U 2 m s, \F{Q ) TV)-, 9 Q I U 2 m si )l 



Q2 _|_ m 2 



^((u 2 
(u 2 - u' 2 ') 2 + 



W\2> 



l u 2 



U 



"|2 /^( U 2- u 2) 



0™/ 



(A47) 
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i(u' 2 < l4 c) K< )i = i(4< \HQ 2 ) Q 3 Kr< )i 

_ F((u 2 -u 2 ') 2 ) , , „ |3 (uf^') s 

u 2 - u 2l *V m » > (A- 18) 



(u 2 — u 2 ) 2 + m. ; 

and for the different parts of the D-term, see Eq. (|A-3|) . we obtain: 

ill I wrijy) I th 1 II II II \ ill I ita(-D) 2 th I // // // \ 

i(u 2 m Si m S3 \V 3 \ \u 2 m Si m S3 ) 1 = 1 {u 2 m Si m S3 \V^ 1 ' |u 2 m Si m S3 )i 

/ \Tr(B)\ II II 
l si 1^31 >2 m Sl 



<W,„ m" 1 ( u 2 m' I y 3 ( 1 B) I u 2 m" ) 1 , ( A- 19) 



1 \ u 2 m si m s 3 I K 31 I u 2 m si m s 3 A 

= j( u' 2 m' si m' S3 |F(Q 2 ) -^-^ Q 2 <x 3 • Q| u' 2 ' m» }:, 



; u' 2 < |V 3 f I u 2 ' r< )x 0^- m u ,P, (A-20) 



"3 s 3 



i(u 2 m si m S3 |V 31 |u 2 m sl m S3 )i 

/II I \fr(D) 5" 1 I // // // \ 

= i( u 2 m si m S3 |V 3 \ I u 2 m si m S3 )i 



= i( u 2 m' si m' S3 \F(Q 2 ) ^2 ^ (°" 3 ' Q) 2 | u 2 m "i m "s >l 

= l(u 2 m' Sl |^K<h 0jf3° 2 , (A-21) 

The matrix elements of V 23 can be evaluated by following the same algorithm 
as above. So we obtain: 

/ /// III \rr{i)\ llll Illl \ £-3/ III llll\ c-3/ III llll\ r 

2 (u a \V 2 \'\u a ) 2 = 6 (u x - u x ) d (u 3 - u 3 ) 5 a m a m, 

x 2 (u 2 a s |y 2 V|u 2 a s ) 2 , (A-22) 

where 

|T/( A ' B ' C ') I J 1 " \ — V r, 2 



/ /// /// | T r(.A,£f,U'J I //// //// \ \^\^2 2 c e £ 

2 (u 2 a s |K 2 3 |u 2 a 5 ) 2 = ^ 2^ 9 a'^'Ys 9 <*s'~is" m '*i TO "i d m%m'^0'w"' 

X 2(u 2 m S2 |y 2 v 3 ; |u 2 m S2 ) 2 , (A-23) 

/ /// /// |ta(-D)i W/ /W \ \ ^ \ A 2 2 r s 

2 {u 2 a s \V 23 |u 2 a s ) 2 = Z^Z^ 9 ^'^' 9a%»T%» °m™m»»o m ™m'J» 



7- 7 »" 



x 2 (u 2 m S2 m S3 \V 23 |u 2 m S2 m S3 } 2 - [A- 24) 

The -B- and C-terms which have been considered in the right side of Eq. 
(|A-23p can be evaluated as: 

/ 111^.111 It/(^)i mi ^.ll'l \ / infill \T?(r\i2\ ^ r\l\ mi I'll \ 

2 (u 2 m S2 \V 2 \ '\u 2 m S2 ) 2 = 2 (u 2 m S2 \F{Q ) Q | u 2 m S2 ) 2 
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F((uf - u' 2 ") 2 ) (uf^') 

I llll ltl\<) i 9 u 2 u 2 VJ m l "m"" 

(u 2 — u 2 Y + m~ m s 2 m s 2 



2 «r< |^W'<>2 = a (« |F(Q' 2 ) Q ,2 |«>2 



(uf - U 2 ") 2 + m 2 , ' 1 Z> 



2 <u' 2 "< \V™ |«> 2 = 2 «r< |F(Q' 2 ) Q' 3 |U2"<'>2 

_ F(« - Ugff |uW/ _ u l,3 («') (A . 27) 

" (uf - u 2 ") 2 + |U2 " 2 1 ' (A 27j 



also the different parts of the D-term which have been considered in the right side 
of Eq. (|A-24|) can be evaluated in the following: 

/ ™W t-fr{D) I th I lilt till lilt \ I in III III \-tAD) 5 th I //////// „«// \ 

2 (u 2 m S2 m S3 | K,^ | u 2 m Sl m S3 ) 2 = 2 (u 2 m S2 m S3 (V^ | u 2 m Sl m S3 ) 2 

= 5m«m™ 2(u 2 m S2 |K 23 |u 2 m si ) 2 , 

(A-28) 



/ /// /// /// ita(^) 2" 1 i //// //// //// \ 

2 ( u 2 m S2 m S3 1^ | u 2 m si m S3 ) 2 

— / ™"' „'» lT/( £) ) 4 th I „//// llll ^llll \ 

- 2 { u 2 m S2 m S3 | V 23 | u 2 m si m S3 ) 2 



/ HI I'! \ f,,\l2, a 2 ' Q / s/2 ,- A',2| "I' "'I: 

2 < u 2 m S2 |F(Q ) Q/2 - m2 Q (0-3 • Q ) | u 2 m S2 ) 2 



2 <u 2 "< l4 B) |uf <>2 05'"5 ) 2 , (A-29) 



<>2 s 2 



2\u 2 m S2 m S3 \V 23 |u 2 m si m S3 , 

— I ■>■>"' ™'" \ T?(n' 2 \ ° 2 ' Q n' 2 ^- c\'\ \ 

- 2 {u 2 m S2 \±< (Q ) ^ n ■ _ 9 Q cr 3 -g |u 2 m S2 ) 2 



Q' 2 + rm 2 

III III i T r(B)i //// nil \ r\( u 2"~ u 2") ' \ on i 

2 (u 2 m S2 \V 2 \ J \u 2 m Sl ) 2 & m i,i m ,J ■ (A-30) 



Appendix B 

Parity and Time Reversal Properties of the Total 4N Wave Function 

In our formulation of 4N bound state we have not yet used the properties of 
total wave function under the parity and time reversal invariance. In this section we 
discuss about these properties. 
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The parity invariance would mean: 

(ua\W) = (ua\PP\&) 

= {-ua\P\¥), P\&) = \¥) 

= (-ua|!?) (B-l) 

The parity invariance or equivalently the symmetry property of the total wave 
function under the exchange of the vector sets |u) to | — u) can be used in the 
numerical solution of the coupled equations (|4-7|) to reduce the dimension of the 4N 
problem. On the other hand this symmetry property can be used explicitly to cut 
down the size of the Yakubovsky components and thus save the time and memory 
when computing the coupled three dimensional integral equations. 

The time reversal invariance might be more interesting. The total wave function 
can be written as: 

( u a \W) = £>| 7 ) ( u 7 \*) = E 3*7 ( u 7 W), (B-2) 

7 7 

where 

(u = (u j\TT\$) 

4 

i=l 
4 

= Y[(-) il ~ ms i~ mt J (-u --y\T\&). (B-3) 

i=i 

By using the time reversal invariance of the total wave function, which is ap- 
plicable for the ground state of 4 He with the total angular momentum J = 0, i.e. 
T\\P) = as well as the parity invariance, Eq. (|B-ip . we can rewrite Eq. f|B-3[) as: 

4 

( u 7 \W ) = Yl (-)( 1 ~ m H- m h) ( u -~f\¥) (B-4) 

i=l 

The time reversal, and the parity invariance, as indicated in Eq. (|B-4p can be also 
used to cut down the size of the total wave function, by reducing the dimension of the 
spin-isospin parts of total wave function, which will be valuable in the calculations 
of 3NFs. 

References 

1) E. Hiyama et al, Phys. Rev. Lett. 85, 270 (2000). 

2) J. Usukura, K. Varga and Y. Suzuki, Phys. Rev. B59, 5652 (1999). 

3) M. Viviani, A. Kievsky and S. Rosati, Phys. Rev. C71, 024006 (2005). 

4) R. B. Viringa et al, Phys. Rev. C62, 014001 (2000). 

5) P. Navratil, J. P. Vary and B. R. Barret, Phys. Rev. C62, 054311 (2000). 

6) N. Barnea, W. Leidemann and G. Orlandini, Phys. Rev. C67, 054003 (2003). 

7) N. W. Schellingerhout, J. J. Schut, and L. P. Kok, Phys. Rev. C46, 1192 (1992). 

8) R. Lazauskas and J. Carbonell, Phys. Rev. C70, 044002 (2004). 



S. Bayegan, M. R. Hadizadeh, and W. Glockle 



H. Kamada and W. Glockle, Nucl. Phys. A548, 205 (1992). 

W. Glockle and H. Kamada, Nucl. Phys. A560, 541 (1993). 

A. Nogga, H. Kamada and W. Glockle, Phys. Rev. Lett. 85, 944 (2000). 

A. Nogga, Ph.D. thesis, Ruhr-Universitat, Buchum (2001). 

H. Kamada et al, Phys. Rev. C 64, 044001 (2001). 

A. Nogga, H. Kamada, W. Glockle and B. R. Barrett, Phys. Rev. C65, 054003 (2002). 

E. Epelbaum et al, Phys. Rev. Lett. 86, 4787 (2001). 

E. Epelbaum et al, Eur. Phys. J. A15 , 543 (2002). 

E. Epelbaum et al, Phys. Rev. C70, 024003 (2004). 

M. R. Hadizadeh and S. Bayegan, Few Body Systems 40, 171 (2007). 

M. R. Hadizadeh and S. Bayegan, Eur. Phys. J. A36, 201 (2008). 

I. Fachruddin, Ch. Elster, and W. Glockle, Phys. Rev. C68, 054003 (2003). 
S. Bayegan, M. Harzchi and M. R. Hadizadeh, in preparation. 

S. Bayegan, M. R. Hadizadeh and M. Harzchi, to appear in Few Body Systems. 
\arXiv:0711403E\ 

S. Bayegan, M. R. Hadizadeh and M. Harzchi, Phys. Rev. C77, 064005 (2008). 
R. Machleidt, Adv. Nucl. Phys. 19, 189 (1989). 

I. Fachruddin, Ch. Elster, and W. Glockle, Phys. Rev. C62, 044002 (2000). 
S. A. Coon, W. Glockle, Phys. Rev. C23 (1981) 1790. 
S. A. Coon et al, Nucl. Phys. A317 (1979) 242. 

M. E. Rose, Elementary Theory of Angular Momentum (Wiley, New York, 1957). 

D. Hiiber, H. Witala, A. Nogga, W. Glockle and H. Kamada, Few Body Systems 22, 107 

(1997). 



